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Abstract—Sangamagrama Madhava (c.1350 - 1425 CE), the 
founder of the Kerala School of Astronomy and Mathematics, 
is credited with the discovery of the power series expansions 
of the trigonometric sine and cosine functions more than two 
centuries before these series were rediscovered in Europe by Sir 
Isaac Newton and Gottfried Wilhelm Leibniz in c.1665 CE. This 
paper examines and analyzes Madhava’s contribution from an 
algorithmic point of view and shows that modern algorithms 
for the computation of the values of sine and cosine functions 
employ essentially the same ideas and methods used by Madhava 
for the computation of the values of these functions. For this 
purpose Madhava’s algorithms have been compared with modern 
computer programmes for the evaluation of these functions 
which are included in standard run time C libraries. The paper 
also observes that Madhava had used Horner’s scheme for the 
evaluation of polynomials, a fact that has not been noticed 
by researchers so far. The descriptions of the power series 
expansions of the sine and cosine functions in Madhava’s own 
words and also their rendering in modern notations are discussed 
in this paper. The sine table computed by Madhava has also been 
reproduced. The paper concludes by stressing the necessity of 
making further investigations into the contributions of the Kerala 
School of Astronomy and Mathematics from an algorithmic point 
of view. 

This paper is also a tribute to the computational genius of 
the Kerala School which flourished in a geographical area very 
proximal to the venue of this International Conference. 


I. INTRODUCTION 


HE trigonometric sine and cosine functions are two ele- 

mentary transcendental functions in mathematics which 
are extensively used in all areas of science and engineering. 
Development of methods for the accurate and fast computation 
of the values of these functions is an area of much interest to 
numerical analysts and computationalists. In this context it is 
highly interesting to explore how this problem was tackled 
in pre-modern times. It would come as a great surprise for 
many that nearly six centuries ago, even before Sir Isaac 
Newton was born, mathematicians in remote villages in the 


Indian state of Kerala had developed computational schemes 
essentially similar to modern algorithms for the computation of 
the values of these functions. Sangamagrama Madhava (c.1350 
- c.1425 CE) was the pioneer of these developments. He 
discovered the power series expansion of the sine and cosine 
functions, developed the Horner’s scheme for the computation 
of polynomials nearly four centuries before the scheme was 
invented in Europe, and successfully applied these results to 
generate a table of values of the sine function accurate to about 
eight places of decimals. 

This paper discusses these developments. After this intro- 
ductory Section I, in Section II of the paper we briefly consider 
the historical context of these ideas. We have have introduced 
the main characters of the Kerala School of Astronomy and 
Mathematics and also the geographical area where they had 
flourished. In Section III we have explained some of the 
relevant terminology used by pre-modern mathematicians of 
India. This includes a discussion of the Katapayddi scheme 
which is a method for the representation of numbers using 
letters of the Sanskrit alphabet. Section IV is an introduction 
to Madhava’s method for the computation of the value of 7. 
In Section V, VI and VII, we present English translations of 
Sanskrit passages containing descriptions of the power series 
expansions of sine and cosine functions and their rendering 
in modern notations. In section VIII we have reproduced 
Madhava’s sine table. Section IX is devoted to an analysis of 
Madhava’s algorithm from a modern perspective. In Section X, 
we compare Madhava’s algorithm for the computation of sine 
and cosine functions with the current standard schemes for 
the computation of these functions. We conclude the paper in 
Section X1 by stressing the necessity for more intensive study 
of these ideas which is likely to lead to better insights into the 
understanding of modern computational schemes. 

In the transliteration of the names and words originally 
given in Sanskrit, we have followed ISO 15919 standard 
which is an extension of International Alphabet of Sanskrit 
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Transliteration (IAST). 


II. KERALA SCHOOL OF ASTRONOMY AND MATHEMATICS 


Even though the contributions of the mathematicians of 
the “classical period” of Indian mathematics are well known, 
the contributions of mathematicians of a later period are 
not so well known (see, for example, [1]). Many historians 
of mathematics were of the belief that generation of new 
mathematical ideas in India generally ceased to exist after 
Bhaskara II. However, research in the past half century has 
revealed the flourishing of a remarkable galaxy of mathemati- 
cians and astronomers in Kerala who had made substantial 
contributions to the development of mathematics. C.M. Whish, 
a Civil Servant with the East India Company, had already 
reported the existence of such a group of personages and 
their accomplishments. (See [2]. This has been reproduced in 
[3].) But his paper was generally ignored by scholars. These 
mathematicians and astronomers, who lived over a period from 
the fourteenth to eighteenth century CE are now collectively 
referred to as forming the Kerala School of Mathematics. 

For a detailed discussion of the mathematical and astronom- 
ical accomplishments of the Kerala School one may refer to 
[1], [4] and [5]. 


A. Prominent members of the Kerala School of astronomy and 
Mathematics 


It is generally agreed that Sangamagrama Madhava who 
flourished during c.1350 - c.1425 CE is the founder of this 
School. No personal details about Madhava has come to light 
(see [6]). The appellation Sangamagrama is surmised to be 
a reference to his place of residence and some historians has 
identified it as modern day Irinjalakuda in Thrissur District 
in Kerala. Only one or two minor works of Madhava have 
survived. But there are copious references and tributes to 
Madhava in the works of other authors about whom accurate 
details are available. It is these references that establish con- 
clusively that Madhava was the founder of the Kerala School. 

Paramesvara (c.1380 - c.1460), a pupil of Sangamagrama 
Madhava, was instrumental in introducing the Drigganita 
system of astronomical computations in Kerala. Damodara, 
another prominent member of the Kerala school, was his 
son and also his pupil. Paramesvara was also a teacher of 
Nilakantha SOmayaji (1444 - 1544), 

One of the most influential works of Nilakantha SOmayaji 
was the comprehensive astronomical treatise Tantrasamgraha 
completed in 1501 (see [7]). He had also composed an 
elaborate commentary on Aryabhatiya called the Aryabhatiya 
Bhasya. In this Bhasya, Nilakantha had discussed infinite 
series expansions of trigonometric functions and problems 
of algebra and spherical geometry. Grahapareeksakrama is 
a manual on making observations in astronomy based on 
instruments of the time. 

Jyésthadéva (c.1500 - c.1600) was another prominent mem- 
ber of the Kerala School. He is most famous for his work titled 
Yuktibhdsa (see [8]). The work was unique in several aspects. 
Firstly, it was composed in Malayalam a regional language 
of the Indian state of Kerala. This was in opposition to the 


then prevailing practice of writing technical works in Sanskrit, 
the language of scholarship. Secondly, it contained elaborate 
proofs of the results expounded in the treatise. Incidentally 
this treatise proves that the idea of proof is not alien to the 
mathematical traditions of India. Yuktibhdasa also contains clear 
statements of the power series expansions of the sine and 
cosine functions and also lucid expositions of their proofs. 

Citrabhanu (c.1550), Sankara Variar (c.1500 - c.1560, au- 
thor of Kriya-kramakari), Acyuta Pisarati (1.1550 1621), 
Putumana Somayaji who authored Karanapadhati and Sankara 
Varman (1774 - 1839, author of Sadratnamdla) were some of 
the subsequent prominent members of the Kerala School of 
Astronomy and Mathematics. 


B. Geographical area of the Kerala School 


As has already remarked, the place where Madhava, the 
founder of the Kerala School, flourished is not known ex- 
actly. However, it is known with certainty that Parameésvara, 
Damodara, Nilakanta Somayaji, Jyésthadéva and Acyuta 
Pisarati all flourished in an area around modern day Tirur, 
Trikkandiyur, Alathiyur and Kudallur in the Malppuram Dis- 
trict in Kerala. These are not far away from Kuttippuram. (See 
Figs. 1.) 


III. SOME TERMINOLOGY 
A. Jya, kdjya, Sara, capa 

Jya, koti-jya (or köjyā) and utkrama-jya (or Sara ) (see 
Fig. 2) are three trigonometric functions introduced by Indian 
astronomers and mathematicians. The earliest known Indian 
treatise containing references to these functions is Siirya Sid- 
hanta. These are functions of arcs of circles and not functions 
of angles. jya and koti-jy are closely related to the modern 
trigonometric functions of sine and cosine. In fact, the origins 
of the modern terms of sine and cosine has been traced back 
to the Sanskrit words jya and k6ti-jya. 

Let ‘arc AB’ denote an arc whose two extremities are A 
and B of a circle with centre O. If a perpendicular BM be 
dropped from B to OA, then (see Fig 2): 

e jyā of arc AB = BM 

e koti-jyā of arc AB = OM 

e utkrama-jyd (or Sara of arc AB) = MA 

If the radius of the circle is R and the length of arc AB is 
s, the angle subtended by arc AB at O measured in radians 
is 0 = $. Then we have (see Fig. 2): 

e jya (arc AB) = Rsin ($). 

e ko-jyd ( arc AB ) = Rcos ($) 

e utkrama-jyd (arc AB)=R (1 — cos (4)) 


B. Katapayddi scheme 


The Katapayadi scheme is a method for representing num- 
bers using letters of the Sanskrit alphabet (see [9]). The 
inventor of this system of numeration is not known. But it had 
been used by Kerala astronomers and mathematicians as early 
as fifth century CE. The scheme was very popular among all 
sections of society in Kerala especially among astronomers 
and mathematicians. The system was virtually unknown in 
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Fig. 1. The geographical area where the Kerala School of Mathematics and Astronomy flourished 


vyasGrdha capa 
(radius) (arc) 
koti-jya utkrama-jya 


(or kojyā) (or Sara) 


Fig. 2. Diagram showing cdpa, jyd , etc. 


North India. The katapayddi scheme is employed in Karantic 
musical traditions to code the serial number of a melakarta 
raga. 

There are several variations of the scheme. In one such 
scheme popular among members of the Kerala School of 
Astronomy and Mathematics, the assignment of digits to letters 
is as in Table 1. 

The following procedure was followed in representing large 
numbers in the katapayddi scheme. 

e Consonants have numerals assigned as per Table I. For 
example, ba is always 3 whereas 5 can be represented by 
either ħa, na or ma or Sa. 

e All stand-alone vowels like a and i are assigned to 0. 

e Incase of a conjunct, consonants attached to a non-vowel 


TABLE I 
MAPPING OF DIGITS TO LETTERS IN katapayddi SCHEME 


ta | tha | da | dha | na | ta | tha | da | dha | na 


pa | pha | ba | bha | ma | - - - - - 


will be valueless. For example, kya is formed by k + y 
+ a. The only consonant standing with a vowel is ya. So 
the corresponding numeral for kya will be 1. 
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e Numbers are traditionally written in increasing place 
values from left to right. The number 386 which denotes 
3 x 100+8 x 10+6 in modern notations would be written 
as 683 in pre-modern Indian traditions. 
Thus Sadratnamdala gives the value of m to 18 decimal 
places in the katapayddi scheme as follows (see [9]): 
bha-dra-mbu-ddhi-si-ddha-ja-nma-ga-ni-ta-Sra- 
ddha-sma-ya-d bhii-pa-gi 
This denotes 423979853562951413 in Indian traditions and 
314159265358979324 in modern notations. 


IV. MADHAVA’S VALUE FOR 7 


To complete the numerical computations of the values of 
the sine and cosine functions, Madhava needed accurate value 
of the mathematical constant 7. The value of m mentioned 
in Sadratnamdla which is correct to 18 decimal places and 
which had been referred to in Section III-B was not available 
to Madhava. Madhava had derived the following series for the 
computation of pi: 


This series is known in European mathematical traditions 
as the Gregory series and its discovery has variously been 
attributed to Gottfried Wilhelm Leibniz (1646 - 1716) and 
James Gregory (1638 1675). The series was known in Kerala 
more than two centuries before its European discovers were 
born (see [10]). Using this series and several correction terms 
Madhava computed the following value for m (see p.119 [11]): 


m = 3.1415926535922. 


V. MADHAVA SERIES FOR SINE AND COSINE FUNCTIONS 


None of Madhava’s works containing any of the series 
expressions attributed to him has survived. These series ex- 
pressions are found in the writings of the followers of Madhava 
in the Kerala school. At many places these authors have clearly 
stated that these are “as told by Madhava”. Thus the enunci- 
ations of the various series found in Tantrasamgraha and its 
commentaries can be safely assumed to be in “Madhava’s own 
words”. The translations of the relevant verses as given in the 
Yuktidipika commentary of Tantrasamgraha (also known as 
Tantrasamgraha-vyakhya ) by Sankara Variar (c.1500 - c.1560 
CE) are reproduced below. These are then rendered in current 
mathematical notations (see [11]). 


VI. MADHAVA’S SINE SERIES 
A. In Madhava’s own words 


Madhava’s sine series is stated in verses 2.440 and 2.441 
in Yukti-dipika commentary (Tantrasamgraha-vyakhya) by 
Sankara Variar. A translation of the verses follows (see p.114 
[11]. 

Multiply the arc by the square of the arc, and take 
the result of repeating that (any number of times). 
Divide (each of the above numerators) by the squares 
of the successive even numbers increased by that 
number and multiplied by the square of the radius. 


Place the arc and the successive results so obtained 
one below the other, and subtract each from the 
one above. These together give the jiva, as collected 
together in the verse beginning with “vidvan” etc. 


B. Rendering in modern notations 
Let r denote the radius of the circle and s the arc-length. 


e The following numerators are formed first: 
S: 8°, S: 8^ S, S- 8f Sf 8",- 


e These are then divided by quantities specified in the verse. 


Ss 
°° 2 + 2)r2? 
s2 s2? 
°° 24 Q)r2 | (42 + 4)r?’ 
s2 s? s2 
°° Fr? FIr (62 46)r2’ 


e Place the arc and the successive results so obtained one 
below the other, and subtract each from the one above to 
get the following expression for jiva: 


os f S 
jiva =s- |s- CERNE 
- 2 2 
S . . 
L (22 +2)r? (42 + 4)r? 
i 8 2 s2 
[$ (22 +2)r2 (42 +4)r2 (62+ 6)r? 


C. Transformation to current notation 


Let 0 be the angle subtended by the arc s at the center of 
the circle. Then s = rô and jiva = rsin. Substituting these 
in the last expression and simplifying we get 

. 6 | g 

EE a 
which is the infinite power series expansion of the sine 
function. 


D. Mādhava’s reformulation for numerical computation 


The last line in the verse “as collected together in the verse 
beginning with vidvān etc.” is a reference to a reformulation 
of the series introduced by Mādhava himself to make it conve- 
nient for easy computations for specified values of the arc and 
the radius (see p.116 - 117 [11]). For such a reformulation, 
Madhava begins by considering an arc of a circle consisting of 
one quarter of a full circle. The measure of a quarter of a circle 
is taken as 5400 minutes (say C minutes). He then develops a 
scheme for the easy computations of the jiva-s of the various 
arcs of such a circle. Let R be the radius of a circle one 
quarter of which measures C. Madhava had already computed 
the value of 7 using his series formula for 7. Using this value 
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of 7, namely 3.1415926535922, the radius R is computed as 
follows: 


R = 2 x 5400/7 
= 3437.74677078493925 
= 3437 arcminutes 
44 arcseconds 
48 sixtieths of an arcsecond 
= 3437 44” 48” 


The expression for jiva can now be put in the following form: 


s’ 8° 


ROF) RF AHA) 
(3) APART 
= see 
C 3! C 5! 
Madhava now pre-computes the values given in Table II. 


With these values, the computational scheme for the evaluation 
of the sine function takes the following form: 


jiva =s 


)?[(2220'39"40'”) 
VICI 5? 47” 
C)?[(16'05"41"") 
)?[(33”0 6”) 
)?(44" I] 


VII. MADHAVA’S COSINE SERIES 


A. In Madhava’s own words 


Madhava’s cosine series is stated in verses 2.442 and 2.443 
in Yukti-dipika commentary (also known as Tantrasamgraha- 
vyākhya) by Sankara Variar. A translation of the verses follows 
(see p.115 [11]). 


Multiply the square of the arc by the unit (ie. 
the radius) and take the result of repeating that 
(any number of times). Divide (each of the above 
numerators) by the square of the successive even 
numbers decreased by that number and multiplied 
by the square of the radius. But the first term is 
(now)(the one which is) divided by twice the radius. 
Place the successive results so obtained one below 
the other and subtract each from the one above. 
These together give the sara as collected together 
in the verse beginning with stena, stri, etc. 


B. Rendering in modern notations 


Let r denote the radius of the circle and s the arc-length. 


e The following numerators are formed first: 


e These are then divided by quantities specified in the verse. 


2 


Ss 
r es 
(22 — 2)r’ 
s? s? 
Dr Ar? 
s2? s2 s2 
are are 66r 


e Place the arc and the successive results so obtained one 
below the other, and subtract each from the one above to 
get the following series for fara: 


8 
sara =r: (2 = Dr? 
f- l s? s? 
(2-22 -r 
[r ; s2 s? s2 
age Uaa 0) 


C. Transformation to current notation 


Let 0 be the angle subtended by the arc s at the centre of 
the circle. Then s = r0 and sara = r(1 — cos0). Substituting 
these in the last expression and simplifying we get 


o o ae 


1- cos0= z- g at 


which gives the power series expansion of the cosine function. 


D. Mādhava’s reformulation for numerical computation 


The last line in the verse “as collected together in the 
verse beginning with stena, stri, etc.” is a reference to a 
reformulation introduced by Madhava himself to make the 
series convenient for easy computations for specified values 
of the arc and the radius (see p.118 [11]). As in the case 
of the sine series, Madhava considers a circle one quarter of 
which measures 5400 minutes (say C minutes) and develops a 
scheme for the easy computations of the sara-s of the various 
arcs of such a circle. Let R be the radius of a circle one quarter 
of which measures C. Then, as in the case of the sine series, 
Madhava gets R = 3437/4448". 

Madhava’s expression for sara corresponding to any arc s 
of a circle of radius R is equivalent to the following: 

8? st 


R?(2? — 2) R- pap = 
A (BG) e] 


Next Mādhava pre-computes the values given in Table III. 
With these values, the computational scheme for the evaluation 


jiva = R. 
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TABLE II 
MADHAVA’S PRE-COMPUTED COEFFICIENTS IN THE POWER SERIES FOR THE SINE FUNCTION 
No. | Expression Value In katapayddi system 
1 R x (n/2)3/3! 2220’ 39” 40" | ni-rvi-ddhd-nga-na-ré-ndra-rung 
2 R x (1/2)°/5! 273’ 57” 47” sa-rvā-rtha-śsi-la-sthi-ro 
3 | Rx (7/2) /7! 16’ 05” 41” ka-vi-Sa-ni-ca-ya 
4 R x (1/2)9/9! 33” 06” tu-nna-ba-la 
5 | Rx (7/2)!+/11! 44" vi-dvān 
TABLE III 
MADHAVA’S PRE-COMPUTED COEFFICIENTS IN THE POWER SERIES FOR THE COSINE FUNCTION 
No. | Expression Value In katapayddi system 
1 R x (7/2)? /2! 424109” 00” u-na-dha-na-krt-bhu-re-va 
2 R x (1/2)*/4! 87203” 05” mi-nā-ngo-na-ra-sim-ha 
3 R x (1/2)°/6! 07143” 24” bha-drā-nga-bha-vyā-sa-na 
4 R x (n/2)8/8! 03°09” 37” su-ga-ndhi-na-ga-nud 
5 | Rx (1/2)!°/10! 05/12!” stri-pi-Su-na 
6 | Rx (7/2)!?/12! 06” ste-na 


of the cosine function takes the following form: 


sara = (s/C)°[(4241'09”00”") 
— (s/C)?[(872'03'"05"”) 
— (s/C)?[(071'43"24"”) 
— (s/C)?[(03'09'"37'”) 
— (s/C)*[(05"12"") 
— (s/C)? 06") II] 


VII. MADHAVA’S SINE TABLE 


The table of values of the sine function prepared by 
Madhava is given in Table IV (see p.121 - 122 [11]). The 
entries in the fourth column of the Table are calculated using 
the entries in the third column of he Table. Let the eight 
digits in the third column of a row corresponding to an angle 
A be dj,d2,...,dg read from left to right. According the 
conventions of the katpayddi scheme, the corresponding value 
in the fourth column is computed as below: 


T 
FD gg |(1000ds + 100dr + 10de + ds) 
10d, +d3 _ 10d) +d 
© 60 3600 


IX. ANALYSIS OF THE ALGORITHMS 
A. These are indeed algorithms! 


In the traditional accounts of the history of Kerala math- 
ematics emphasis is placed on the discovery of the power 
series expansions of the sine and cosine functions. The fact 
that the Kerala School had produced acceptable and sound 
proofs of these series expansions is also strongly emphasised. 
The methods and techniques used in the proofs have generated 
much interest and more controversy among the historians 


of mathematics, some arguing that these are the beginnings 
of calculus and Madhava and his Kerala School should be 
accepted as the inventors of calculus while many others 
propounding that though the proofs are basically correct, 
the methods do not conclusively establish the invention of 
calculus. But, one fact that emerges from a close reading of the 
relevant passages is that the Kerala inventors of these series 
were more interested in the procedures that could be developed 
based on these series. For them these series provided nothing 
but an algorithm for the computation of the sine and cosine 
functions. This algorithmic aspects are all the more clear in 
the way the series expansions were formulated. It is in fact 
given as a step by step procedure for easy computations the 
function values. Moreover in Jyésthadéva’s Yuktibhasa, the 
author has used precisely the term kriya-krama which clearly 
and unambiguously translates into procedure or, in modern 
terminology, an algorithm. Madhava and other members of the 
Kerala School were fundamentally computationalists. Their 
goal was to compute the values of sine and cosine functions. 
For this they developed the necessary mathematics and then 
created the algorithm. The algorithm was then actually applied 
to create a table of values of the sine function correct to 
about eight decimal places. The procedures outlined above are 
Madhava’s algorithms for the computations of sine and cosine 
functions. 


B. Use of a aolynomial approximation 


Madhava’s algorithms for the computation of sine cosine 
functions use polynomial approximations. The algorithm for 
the computation of the sine function makes use of an 11 
th degree polynomial where as the algorithm for the cosine 
function makes use of a 12 th degree polynomial. The orders of 
the polynomials were decided by the requirements of accuracy. 
They were targeting an accuracy to a third (that is, one- 
sixtieth of a second) which would suffice for their astronomical 
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TABLE IV 
MADHAVA’S SINE TABLE 
Angle A Madhava’s numbers Madhava’s numbers | Value of sin A Modern value 
in degrees | for specifying sin A for specifying sin A | derived from of sin A 
in ISO 15919 in Arabic numerals Madhava’s table 
transliteration scheme 

03.75 Sresth6é nama varisthaném 22 05 4220 0.065403 14 0.065403 13 
07.50 himadrirvédabhavana: 85 24 8440 0.13052623 0.13052619 
11.25 tapan6 bhanu siiktajno 61 04 0760 0.19509032 0.19509032 
15.00 maddhyamam viddhi dohanam | 51 54 9880 0.25881900 0.25881905 
18.75 dhigdjy6 nasanam kastam 93 10 5011 0.32143947 0.32143947 
22.50 channabhogasayambika 70 43 5131 0.38268340 0.38268343 
26.25 mrgaharo narésoyam 53 82 0251 0.44228865 0.44228869 
30.00 virð ranajayotsuka: 42 25 8171 0.49999998 0.50000000 
33.75 müla: visuddha: ndlasya 3 45 9091 0.55557022 0.55557023 
37.50 gdnesu viralā nara: 30 64 2902 0.60876139 0.60876143 
41.25 asuddhigupta corasri: 05 93 6622 0.65934580 0.65934582 
45.00 Samkukarno nagésvara: 51 15 0342 0.70710681 0.707 10678 
48.75 taniijo garbhaj6 mitram 60 83 4852 0.75183985 0.75183981 
52.50 $rimanatra sukhi sakhé 25 02 7272 0.79335331 0.79335334 
56.25 §aSi ratrou himaharou 55 22 8582 0.83146960 0.83146961 
60.00 végajna: pathi sindhura: 43 01 7792 0.86602543 0.86602540 
63.25 chaya layō gajé nil 71 31 3803 0.89687275 0.89687274 
67.50 nirmald nasti salkulé 05 30 6713 0.92387954 0.92387953 
71.25 ratrou darpanamabhramgam 22 81 5523 0.94693016 0.94693013 
75.00 ndgastumga nakho bali 03 63 0233 0.96592581 0.96592583 
78.75 dhiré yuva kathalola: 92 14 1733 0.98078527 0.98078528 
82.50 pūjyð narijanairbhaga: 11 02 8043 0.99144487 0.99144486 
86.25 kanyagaré nagavalli 11 32 0343 0.99785895 0.99785892 
90.00 devo visvasthali bhrgu: 84 44 7343 0.99999997 1.00000000 


applications. The values computed by Madhava could also 
be obtained by other specialized methods. But Madhava did 
seek and get a general method speaks volumes about his 
computational genius. 


C. Pre-computation of the coefficients 


The observation that Madhava was a computationalist par 
excellence would be reinforced if we note that he did not stop 
with inventing the power series expansions. He also developed 
an efficient computational scheme. One of the ingredients 
of this scheme was the pre-computation of the coefficients 
appearing in the polynomial approximations. The coefficients 
of both sine and cosine series were pre-computed. These pre- 
computed coefficients are given in Table II and Table II. Of 
course these are not in the modern decimal notations. They 
are in the methods and traditions then prevalent in Kerala. 


D. Use of Horner’s scheme 


A final observation is that Madhava had applied what is 
now known as the Horner’s scheme for the computation 
of polynomials. The scheme is now attributed to William 


George Horner (1786 1837) who was a British mathematician 
and schoolmaster (see p.271 [12]). In modern notations, the 
scheme can be formulated as follows. 


Given the polynomial 


n 
p(z) = 5 aix = ao + aya + a22? + 032° +- tana”, 
i=0 


where ao, . . . , @ are real numbers let it be required to evaluate 
the polynomial at a specific value of x. To do this Horner 
suggested that the polynomial be expressed in the form 


p(t) = ao + z(a + z(a +--+ + 2(dn-1 + Gnx)---)). 


This form of the polynomial expression suggests the following 
algorithm for the computation of its value for a given value 
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of x: 


bn = an 


bn—1 = ûn-1 T bns 


bn—2 = An—2 T bp—12 


by = ai T box 


bo = ag T biz. 


bo is the required value of the polynomial. 

Horner’s scheme requires only n multiplications and n 
additions to compute the value of a n th degree polynomial. 

A quick glance at the expressions for sine and cosine using 
the pre-computed coefficients readily shows that Mādhava had 
actually implemented Horner’s scheme in his algorithms. 

Even though the algorithm is named after William George 
Horner, who described it in 1819, the method was already 
known to Isaac Newton in 1669, the Chinese mathematician 
Qin Jiushao in the 13th century, and even earlier to the Persian 
Muslim mathematician Sharaf al-Dn al-Ts in the 12 th century. 
The use of the scheme has also been traced to a Chinese work 
of the 3 rd century. However Madhava’s was the first conscious 
and deliberate application of the scheme in a computational 
algorithm with the intention of reducing the complexity of 
numerical procedures. 


E. Simultaneous computation of sine and cosine 


In many modern implementations of routines for the calcu- 
lations of the sine and cosine functions, there would be one 
routine for the simultaneous computation of sine and cosine. 
The argument in favour of such an approach is that in most 
engineering and other applications, whenever sine or cosine 
is required, the other would also be required. So a common 
algorithm which returns both values simultaneously would be 
more time efficient and economical especially if the main 
component of a certain computational process is the evaluation 
of sine and cosine functions. It would appear that Madhava 
had anticipated such a scenario. This is evidenced by the 
description of one common procedure for the evaluation of 
sine and cosine functions in Yuktibhdsa. We quote a translation 
of this common procedure below (see p.102 [13]): 

This being the case, the following is the procedure 
(kriya-krama) to be adopted. The required is the 
first result. When this squared, halved and divided 
by the radius, the second result is got. Keep this 
second result separately. Now multiply the second 
result also by the arc and divide by 3 and also by 
radius. Place the result got below the first result. 
Then multiply this also by the ar and divide by four 
and the radius. Keep the result below the second 
result. in this manner, derive successive results by 
multiplying the previous result by the arc and di- 
viding by corresponding successive numbers 1,2,3, 
etc. and by radius. Now place below the first result 
the odd results, viz. the third, the fifth, etc. and 
place below the second result the even results, viz, 
the fourth, the sixth, etc. Then subtract successively 


from the bottom result from the one above it, the 
remainder from the one still above it. Ultimately in 
the first column the resultant first result will be left 
and in the second column the resultant second result 
will be left. These will be the required Rsine (jya) 
and Rversine (Sara). 


X. COMPARISON WITH A MODERN ALGORITHM 


A comparison of Madhava’s schme with a modern scheme 
for the computation of sine and cosine values is worth consid- 
ering. One such algorithm openly available in the internet is 
the programme sincos.c included in many Linux distributions. 
For discussion we have taken the programme included in the 
Open64 Compiler developed by Computer Architecture and 
Parallel Systems Laboratory in University of Delaware (see 
[14]). One can also compare Madhava’s scheme with other 
codes for the computation of these functions, for example, the 
codes given in the Netlib Repository maintained by AT&T 
Bell Laboratories, the University of Tennessee and Oak Ridge 
National Laboratory (see [15], [16]). These programs are not 
using the polynomials used by Madhava. Instead they are 
using the minimax polynomial computed using the Remez 
algorithm to improve the accuracy of computations. But the 
fact remains that these modern algorithms are indeed using 
polynomial approximations for computations of the sine and 
cosine functions. Also, they are using pre-computed coeffi- 
cients and Horner’s scheme for the evaluation of polynomials. 
Since the polynomials used are the minimax polynomials, the 
pre-computed coefficients given in these computer programs 
are different from the pre-computed coefficients employed by 
Madhava. 


A. Coefficients for evaluation of sine 


The following program segment specifies the pre-computed 
coefficients in the polynomial approximation for the sine 
function. The values are given in the IEEE:754 floating point 
format. 


00134 

00135 /* coefficients for polynomial 
approximation of sin on +/- pi/4 x/ 

00136 

00137 static const du S[] = 

00138 

00139 D(0x3ff00000, O0x00000000), 

00140 D(0xbfc55555, 0x55555548), 

00141 D(0x3£811111, Ox1110f7d0), 

00142 D(O0xbf2a01la0, 0x19bfdf03), 

00143 D(0x3ec71de3, 0x567d4896), 

00144 D(O0xbe5ae5e5, 0xa9291691), 

00145 D(0x3de5d8fd, Oxlfcf0ecl), 

00146 ; 

00147 
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B. Coefficiets for the evaluation of cosine 


The following program segment specifies the pre-computed 
coefficients in the polynomial approximation for the cosine 
function. The values are also given in the IEEE:754 floating 
point format. 


00148 /* coefficients for polynomial 
approximation of cos on +/- pi/4 */ 
00149 


00150 static const du C[] = 
00151 

00152 D(0x3ff00000, 0x00000000), 
00153 D(Oxbfdfffff, Oxffffff96), 
00154 D(0x3fa55555, 0x5554f0ab), 
00155 D(Oxbf56cl6c, O0x1l640aaca), 
00156 D(0Ox3efa019f, Ox81lcb6ald), 
00157 D(O0xbe927df4, 0x609cb202), 
00158 D(0x3e21b8b9, 0x947ab5c8), 
00159 ; 


C. Polynomial approximations using Horner’s scheme 


The following lines in the code describes the computations 
of the polynomial approximations for the sine and cosine 
functions simultaneously. 


00329 xsq = X*x; 


00330 

00331 cospoly = (((((C[6].d*xsq + 
C[5].d)*xsq + 

00332 C[4].d)*xsq + C[3].d)*xsq + 

00333 C[2].d)*xsq + C[1].d)*xsq + 
C[0].d; 

00334 

00335 sinpoly = (((((S[6].d*xsq + 
S[5].d)*xsq + 

00336 S[4].d)*xsq + S[3].d)*xsq + 

00337 S[2].d)*xsq + S[1].d)*(xsq*xx) + x; 


D. Comparison 


It is true that this programme having more than 500 lines 
of code has made use of several other ideas as well. But the 
critical components continues to be the following which are 
essentially the ideas enshrined in Madhava’s computational 
scheme developed more than six centuries ago. 


e Use of an approximating polynomial. 

e Pre-computation of coefficients. 

e Use of Horner’s scheme for the evaluation of polynomi- 
als. 

e Simultaneous computation of sine and cosine functions. 


XI. CONCLUSION 


Sangamagrama Madhava, the founder of Kerala School of 
Astronomy, was a great computationalist and an algorithmist. 
When confronted with a numerical problem he created the 
necessary mathematical tools with rigorous proofs, devel- 
oped efficient computational schemes, and implemented those 


schemes. Mathematicians and astronomers have studied the 
contributions of the Kerala School in great detail. But it 
appears that his works have not been thoroughly scrutinised 
by algorithmists and computationalists. There is a need of a 
more intensive study of the methods of Kerala School from 
an algorithmic point of view. This will help us understand our 
own past much better. More importantly it may help us develop 
new algorithmic paradigms which may lead to the creation of 
newer computational tools. 
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